Strong terahertz response in bilayer graphene nanoribbons 
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We reveal that there exists a class of graphene structures (a sub-class of bilayer graphene nanorib- 
bons) which has unusually strong optical response in the terahertz (THz) and far infrared (FIR) 
regime. The peak conductance of terahertz/FIR active bilayer ribbons is around two orders of mag- 
nitude higher than the universal conductance of /Ah observed in graphene sheets. The criterion for 
the terahertz/FIR active sub-class is a bilayer graphene nanoribbon with one-dimensional massless 
Dirac Fermion energy dispersion near the F point. Our results overcome a significant obstacle that 
hinders potential application of graphene in electronics and photonics. 
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In recent years, eraphene has attracted a great deal 
of interest[ii[3,lll3. New physics have been predicted 
and observed, such as electron-hole symmetry and half- 
integer quantum Hall effect [1, 0], finite conductivity at 
zero charge-carrier concentrationj2|, and the strong sup- 
pression of weak localization!^ Qj^. Bilayer graphene 
(BLG) has also attracted considerable attention recently, 
with seminal experimental and theoretical work being 
carried out [l|, Q . Bilayer calculations use interlayer cou- 
pling constants based on the Slonczewksi-Weiss-Mclure 
model [ol, |Tc||. By further confining the electrons in 
the graphene plane, one can obtain one-dimensional 
structures which we refer to as graphene nanoribbons 
(GNRs)[lH. It has been shown that GNRs with zigzag 
edges can have finite magnetization with either ferro- 
magnetic order or antiferromagnetic order 12, T^l- These 
properties promise building blocks for technological ap- 
plications in molecular electronic and optoelectronic de- 
vices. 



The optical properties of graphene systems is a topic 
of considerable interest. In particular, the minimal 
conductivity of single layer graphene (SLG) within the 
Dirac regime is a much celebrated result which was 
calculated theoretically long before graphene's fabrica- 
tion in 2003 [ii]. The optical conductivity of SLG 
outside the low ener gy D irac regime has been calcu- 
lated theoretically [islTlill . The conductance of bilayer 
graphene has also been calculated [l3|, as has the con- 
ductance of various GNRs [l^. All of this research 
has shown that the optical responses of graphene and 
graphene nanoribbons are extremely weak. In the EM 
frequency band from terahertz to visible, the absorp- 
tion coefficient for these systems is generally less than 
3%. 1^ 2^ 21 1 . While the lateral confinement and the 
edge states can lead to the ferromagnetic and antiferro- 
magnetic order in GNRs [l^ . the optical response of all 
ribbons remains very weak. There are two fundamental 
reasons for this: (1) the density of states vanishes near 
the Fermi energy, and (2) the interband transition am- 



plitude is small. Because of this, potential application of 
graphene structures in optoelectronics and photonics is 
severely limited. To date, these obstacles have remained. 

However we are now able to demonstrate that there 
exists a sub-class of bilayer graphene nanoribbons 
(BLGNRs) which have an unusually strong optical con- 
ductance in the terahertz (THz) to far infrared (FIR) 
regime. The height of the conductance peak is close to 
two orders of magnitude greater than the universal con- 
ductance of graphene sheets. We found that this sub- 
class of BLGNRs can be either armchair or chiral, but 
their energy dispersion near the F point must be that 
of a one-dimensional massless Dirac Fermion. This sub- 
class of graphene structures are the first systems to show 
such a strong optical response in the absence of any ex- 
ternal field in the important frequency band of THz and 
FIR. 

We first construct a single layer GNR (SLGNR) follow- 
ing the convention of Ezawa [22], i.e., a ribbon is specified 
by two indices p and q. We begin by placing m = p + q 
hexagons next to each other with flat edges touching. On 
top of this layer, we place an identical layer, offset by q 
hexagons. Continuing this in both directions we can con- 
struct a GNR of arbitrary chirality. Within this model 
q = corresponds to a zig-zag (ZZ) edged ribbon, and 
q = I corresponds to an armchair (AC) edged ribbon. 
q > 1 corresponds to arbitrary chirality, defined by an 
angle 9 — tan^^(-\/3/(2g + 1)). The number of atoms 
in the unit cell is given hy Nu = iq + 2p + 2. The sec- 
ond layer is now constructed by assuming the standard 
A — B 'Bernal' stacking, along the vertical C-C vector. 
The net effect is simply shifting the entire second layer 
up (or down) by an amount C — C = 1.42 A. A typical 
BLGNR is shown in Figure 1. 

The intralayer coupling is calculated using the tight 
binding formalism where t « 3eV is the nearest neigh- 
bour hopping integral. The edge effects of GNRs are 
incorporated into the tight binding formalism by setting 
the overlap integral to zero for hopping between edge 
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FIG. 1: The construction of a bilayer graphene nanoribbon. 
The second layer is identical to the first, but shifted along a 
C-C vector. 

sites and their neighbours which are off the edge of the 
ribbon. In this way, the edge states are incorporated into 
the full electronic properties. The interlayer coupling is 
restricted to the dominant coupling term 7 = 0.36eV 
which occurs only between A and B sites which sit di- 
rectly one on top of the other. As an example the Hamil- 
tonian matrix for a (2,1) BLGNR is given by 

^(2,1) =1^* H- , I 
\ -"inter -"intra / 

Where the intralayer terms are given by 
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FIG. 2: The bilayer graphene nanoribbon electronic disper- 
sion curves differ from the single layer ones as each single 
layer subband has a complementary subband which generally 
differs in energy by < 7. The zigzag case is given in (a), and 
the Dirac armchair case in (b). The inset to (b) shows the 
two sets of symmetric bands near the Dirac point. The arrow 
represents transitions between non-symmetric bands which is 
the dominant transition at low energies that leads to the un- 
usually high optical conductance observed in figure 3(b). 
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The band structures for ZZ-BLGNRs are shown in Fig- 
ure 2(a). They differ from the single layer case as the in- 
terlayer coupling produces a second subband correspond- 
ing to each single layer subband offset by an amount 



7 for most of the Brillouin Zone. Near the Dirac-like 
points however, the intersubband gap goes to zero for 
the low energy bands which approach zero energy, and 
the higher energy bands all approach a degenerate point 
slightly away from the Dirac-like points. The edge states 
cause the plateau of the zigzag dispersion at low ener- 
gies, causing an extended region of zero band-gap near 
the zone boundary. 

The band structures of AC-BLGNRs shown in figure 
2(b) behave much like 2D bilayer graphene. The linear 
dispersion at the Dirac points becomes curved, and a 
second subband appears. Each subband is thus paired. 
But unlike the ZZ-BLGNR case, the intersubband gap 
does not go to zero near the Dirac points, but remains at 
least 7/2 separated from it's pair. As the width increases, 
this gap approaches zero as the band pairs converge. This 
is strictly an edge effect, and in the 2D limit a band gap 
of 7 re-emerges. 

We shall show below that the oscillator strength of the 
interband transition is strongly dependent on the prop- 
erties of the energy dispersion at the zero-gap point. We 
shall show below that for ZZ-BLGNR, where the zero gap 
position is at the K point and the two low energy disper- 
sions are close to parabolic, the oscillator strength at low 
energy is very small. On the other hand, for AC-BLGNR, 
the zero gap position is at the T point. Furthermore, the 
two low energy dispersions are very close to linear (or a 
one-dimensional massless Dirac Fermion). In this case 
the interband transition between the bands (transitions 
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FIG. 3: The optical conductivity for the ZZ-BLGNR (a) 
and the AC-BLGNR (b). The low energy activity in the bi- 
layer ribbons is particularly significant, especially in the Dirac 
AC-BLGNR where the optical conductivity is approximately 
80ao. 

between the non-symmetric bands shown by arrows in 
the inset of figure 2(b)) is extremely strong. 

The optical conductivity is calculated using the Kubo 
formula given by 

= dte-*([J(i),J(0)]) (6) 

Where J, the current operator is given by dH/dy. We de- 
fine the dimensionless photon frequency O = hw/t, and 
normalize our results to the single layer universal con- 
ductivity given by cto = e'^/Ah. We determine the depen- 
dence of the optical conductance on the ribbon widths 
and chiralities. In figure 3(a) we show that the opti- 
cal conductivity for the ZZ-BLGNR exhibits a spike cen- 
tered on zero energy. This spike occurs here because 
both low energy subbands approach zero energy at the 
Dirac-like points. In the SLGNR case, the velocity opera- 
tor approaches a constant, which makes inter-symmetric- 
subband transitions forbidden. This is no longer the case 
in bilayer ribbons, and there is also now the possibility 
of low energy inter-non-symmetric-subband transitions. 

Over the full energy spectrum, we see that some of 
the resonant peaks in the single layer optical conductiv- 
ity spectrum have split into three peaks. This will not 
generally be the case. Most peaks will split into two as 
will be seen in the armchair case. However, near $7 = 1, 
the subbands create a linear Dirac-like band structure 
with features similar to the Dirac point in the proceed- 
ing armchair well as those observed in 2D bilayer 
graphene. This means that there are three possible en- 
ergy transitions with high density of states. The central, 
primary peak corresponds to the original SLGNR peak, 
and the two secondary peaks, one below, and and one 



FIG. 4: The width dependence of the energy of the low-energy 
peak for BLGNRs with q = [1,4]. The strongest peaks occur 
in the lowest energy gaps in the Dirac armchair ribbons. The 
inset shows the width dependence of the band gap for Dirac 
BLGNRs with strong low energy optical response. This gap 
eventually disappears, but in the 2D limit with no edges, it 
re-emerges at 7. 

above the original by an amount Q ~ correspond to 
the new curved subbands which don't quite touch the 
degenerate point from the single layer case. 

The optical conductance of armchair BLGNRs is 
shown in figure 3(b), In the AC case it peaks sharply 
at 7 /2 and trails off because of the curvature of the band 
structure. In the armchair peak is still observed 

as J7 ^ 0, but the peak at 7/2 is about 2 orders of mag- 
nitude stronger. This peak corresponds to vertical tran- 
sitions between non-symmetric subbands which are far 
more probable than those between symmetric ones (see 
fig. 5 and subsequent discussion). This single low energy 
peak is larger than every other peak across the spectrum. 
As the width of the ribbons increases, the strength of this 
peak decreases and the width of it broadens. In the case 
of an infinitely wide ribbon, the strong optical response 
is lost, settling at ~ 4(To at = 7. This reflects the 
peculiar edge dependence of the bilayer ribbon system. 

At higher energies, the single peaks observed in AC- 
SLGNRs generally have split into two, and are separated 
by an amount w 27. This corresponds to two sets of sym- 
metric transitions, the non-symmetric transitions being 
largely suppressed. 

For q > I BLGNRs, the band gap between the two 
lowest energy symmetric bands varies from zero to leV. 
Similarly the non-Dirac AC-BLGNRs (ie {p + q)/3^ I), 
have varying band-gaps for the lowest energy subbands. 

For a given type of BLGNR, the conductance peak 
position can be tuned with the ribbon width. Figure 4 
shows the width dependence of the peak position in the 
THz/FIR regime. The peak position oscillates with the 
ribbon width. The amplitude of the oscillation is of the 
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FIG. 5: The p dependence of the low energy peaks for a 
(p,l) BLGNR. In Dirac ribbons the symmetric transitions 
are suppressed and the non-symmetric ones dominate. For 
non-Dirac ribbons the opposite is the case. 



same order of magnitude as the average peak position, 
indicating a large range for tuning the resonance peak. 
The period of the oscillation increases with the chirality 
(q) . The inset of figure 4 shows the width dependence of 
the energy gap for Dirac BLGNRs. This gap decreases as 
width increases making the location of the optical peak 
strongly width dependent. In the limit of infinite width 
this edge determined gap approaches zero, and the 2D 
bilayer band structure sees the emergence of a gap of 7. 
The optical response in this case however is not nearly 
so large, being « Actq 

Figure 5 shows the width dependence of the magni- 
tude of the low energy peak for the {p, 1) BLGNR. For 
(p -|- l)/3 G /, the non-symmetric matrix element dom- 
inates, causing the single low energy peak. The low- 
frequency peak conductance for this class of BLGNRs 
is unusually strong having a value of approximately 
SOctq, much stronger than the universal conductance of 
graphene sheets. 21] When the Dirac condition is not met 
however (ie (p+l)/3 ^ /), the symmetric matrix elements 
dominate, and the non-symmetric matrix elements are 
greatly suppressed, leading to a much weaker response 
to the low energy spectrum. 

The width dependence of the strength of the low en- 
ergy peak for ZZ-BLGNRs is given in Figure 6. For very 
narrow ribbons with p < 6 the peak quickly decreases 
in magnitude at a decreasing rate. For p > 6 however, 
the peak magnitude increases steadily reflecting the low 
energy subband shape. As the width increases, the low 
energy subbands remain lower, which increases the DOS, 
allowing more transitions between subbands. For very 
narrow width ZZ-BLGNRs however, the curvature in the 
subbands is so high that the velocity operator allows 
strong coupling between the subbands, which makes the 



FIG. 6: The width dependence of the magnitude of the low 
energy peak for the ZZ-BLGNRs decreases quickly with in- 
creasing width, and then increases steadily for p > 6. 



low energy magnitude very strong. 

In summary, we have shown that the interplay of rib- 
bon's chirality and the inter-ribbon coupling can lead to 
significant enhancement in optical response. We have 
identified a sub-class of BLGNRs where the inter-ribbon 
coupling causes a finite band gap in the energy minimum 
(maximum) and induces strong inter-subband transi- 
tions. The distinct feature of this sub-class of BLGNRs is 
that they have a one dimensional massless Dirac Fermion 
dispersion near the F point. The peak conductance of 
this class of BLGNRs has a very large value of as much 
as 80(To, making them a class of materials for unique ap- 
plications in optoelectronics. The simple picture behind 
this phenomenon is that the density of states for the ID 
massless Dirac Fermions remains finite at zero energy, 
whereas that for the 2D massless Dirac Fermions in a 
graphene sheet vanishes. The peculiar role of the edge 
states has also been shown to contribute as the 2D limit 
of infinite width shows markedly decreased response in 
this frequency regime. 

These results reported here open a gateway to the cre- 
ation of graphene-based low energy photon devices. The 
ribbon width and chirality selection for various applica- 
tions is crucial, as the optical responses of various ribbons 
change dramatically when these properties are varied. 
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